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Abstract. Let D be an open disk of radius < 1 in C, and let 
(en) be a sequence of ±1. We prove that for every analytic func- 
tion / : _D — )■ C without zeros in D, there exists a unique sequence 
(a„) of complex numbers such that /(z) = /(O) nJ^i(l + ^n^:")"" 
for every z € D. From this representation we obtain a numerical 
method for calculating products of the form JJ^ pnmo fi^/p) P™" 
vided /(O) = 1 and /'(O) = 0; our method generalizes a well known 
method of Pieter Moree. We illustrate this method on a constant of 
Ramanujan Tr"^/^ Hp prime Vp'^ - pin (p/ (p - 1)). From the prop- 
erties of the exponents a„, we obtain a proof of the following con- 
gruences, which have been the subject of several recent publications 
motivated by some questions of Arnold: for every n x n integral 
matrix A, every prime number p, and every positive integer k we 
have tr AP = tr Ap (mod p'^) . 
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1. Introduction 



Many constants in number theory appear in the form Ylpfi^/p)^ 
where the product is taken over all (sufficiently large) prime numbers 
and / is a function analytic in a neighborhood of and such that 
/(O) = 1 and /'(O) = 0. The results of this paper arose from our 
attempt to compute some constants of this type with high accuracy. 
From this perspective, our work should be considered as a generaliza- 
tion of the technique of Pieter Moree |T2], who shows how to compute 
such products to high accuracy when / is a rational function satisfying 
certain additional properties. 

To achieve our goal we prove the following result: 
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Theorem. Let f{z) be an analytic function without zeros in the disk 
1^1 < i? < 1 and let (e„) be a sequence of ±1. There exists unique 
sequence (an) such that the product f{0)ll'^^i{l + 6^2;")"" converges to 
f uniformly on compact subsets of the disk \z\ < R. 

(see Theorem 12.91 where a simple way to compute the exponents a„ 
is stated). As an illustration we mention the following formulas (see 
Examples [SHI Ens]): 



n( 

n=l 



n( 

n=l 



</.(n) 



valid for \z\ < 1, where /i, (j) are the Mobius function and Euler func- 
tion respectively. After this work was completed, Pieter Moree kindly 
pointed to us an interesting paper by G. Dahlquist |6], where a product 
decomposition as in Theorem above is obtained, and it is used to inves- 
tigate analytic continuation of certain Euler products. As a matter of 
fact, the idea to use such product decompositions for various functions 
goes back to works of Mobius [13] and later Landau (see, for example, 
[To]). Dahlquist 's discussion of his result is somewhat brief, therefore 
we provide a detailed proof including a careful analysis of the expo- 
nents a„. We hope that our paper will make this result more widely 
known, which it fully deserves because of its many applications. In 
addition to being instrumental for the numerical method described in 
Section |H it naturally leads to interesting arithmetic and combinato- 
rial results by studying the relations between the exponents and the 
Taylor coefficients of / and some related functions. In Section 3, we 
will show an example of such a result. Namely, we will obtain a short 
proof of the following theorem: if A is an integral matrix, p is a prime 
number, and n is a positive integer then the traces of A^" and A^" 
are congruent modulo p". This result has been conjectured at the be- 
ginning of this century by V. I. Arnold, who considered it as an analog 
of the classical Euler theorem. Arnold's conjecture has been the sub- 
ject of several recent publications (among others, see [2], [H], [TS] , 
[19j . |llj). where several different proofs can be found. Let us mention 
yet another application. The products as in Theorem above have been 
considered as formal identities in the theory of g-series, mainly in the 
case when the exponents a„ are integers (see [1]). Finally, the product 
decomposition has been used in a recent work [Hj, where the authors 
have proved that the regularized product of all prime numbers is equal 
to 47r2. 
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Returning to our original goal to get an accurate approximation of 
the product Yl^ f{l/p), the strategy now is the same as the one devel- 
oped in [12]. Namely, using the decomposition f{z) = 11^^(1 — z"^)""" 
we get that 

OD oo OO 

n /( = n 11(1 - p-^r- = n 11(1 - p-t^- = n ^^^^ ^)""" 

p>t p>t n=2 n=2 p>t n=2 

where t is sufficiently large and ((t,s) = Y[p>t(^ — is a partial 

zeta-function. The key observation behind our method is that the prod- 
uct on the right converges rapidly. As an illustration, in Examples 14.21 
and 14.31 we compute ffist 50 decimal digits of two constants appearing 
in analytic number theory. To the best of our knowledge, this has not 
been done before. After writing a simple code, the computations using 
PARI/GP and an ordinary laptop take only several seconds. 

Acknowledgments. We thank Maxim Korolev for providing us with 
information related to Example 14.31 B. Petrenko thanks Pieter Moree 
for very useful discussions of his work [T2] during Petrenko's visits to 
MPIM in July-August of 2009 and in August 2010. 



2. Product decomposition 

For complex numbers a and z such that \z\ < 1 we write (1 + z)" for 
the binomial series 1 -|- Yl'^=i O'^"- Then (1 + z)°' = e"^"^'^-'^"'"^) , where 
Log(l + z) = is principal branch of logarithm. 

Recall that a sequence (/„) of functions on a topological space X 
converges compactly to a function / on X if it converges uniformly 
to / on every compact subset of X. We will need the following well 
known consequence of the Residue Theorem: 

Lemma 2.1. Let {fn) be a sequence of analytic functions in a domain 
U , none of which assumes the value in U . Suppose that this sequence 
converges compactly on U to a function f . Then either f = or f does 
not assume the value on U. 

Proof. The function / is analytic in U. Suppose that / is not identically 
0. Let zq G U. Then there is e > such that f{z) 7^ for all z ^ zq in 
the disk |z — zqI < It follows that /^//„ converges uniformly to /'// 
on the circle 7 with center zq and radius e. Thus f'{z)dz/ f{z) = 
limn^oo fn{z)dz/ fn{z) = 0, as J^fl^{z)dz/fn{z) = (since /;//„ is 
analytic in U). It follows that f{zQ) 7^ 0, as otherwise the function 
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/'// would have a simple pole at zq and the integral f'{z)dz/ f{z) 
would not vanish. □ 

Proposition 2.2. Consider the product + e„2;")"", where = 

±1, \z\ < 1, and G C. Let < R < 1. The following conditions are 
equivalent. 

(1) The product converges compactly on the circle \z\ < R. 

(2) The series «nLog(l + e„2;") converges compactly on the 
circle \z\ < R. 

(3) The series Yl^=i ^nnanz""^^ / {1 + enZ"") converges for \z\ < R. 

(4) The series Yl'^=i (^nz"' converges for \z\ < R. 

Moreover, if one of the above equivalent conditions holds, then the se- 
ries in (2), (3), and (4) converge absolutely. 

Proof. Assume (1). Then the product defines an analytic function f{z) 
on \z\ < R which does not vanish at any point by Lemma [2.11 . Thus 
h{z) = log/(z) exists and is analytic on \z\ < R (we take here the 
logarithm satisfying log/(0) = 0). Let hn{z) = X^Li afcLog(l + ekz''). 
Let r < R. Then |/(^)| > B for some -B > and all z such that \z\ < r. 
For any e > there is n such that \f{z) — e'^'^^^-*! < Be for all N > n. It 
follows that |1 -e'''v(^)~''(^)| < e for all z such that \z\ < r. For e < 1/2, 
this implies that Log(l — e) < ^{h^{z) — h{z)) < Log(l + e) and 
'^{h]S!{z) — h{z)) G |J^gg(— arccos(l — e) + 2/c7r, arccos(l — e) + 2A;7r). 
Since '^{h^{z) — h{z)) is continuous on the connected set \z\ < R 
and it vanishes at 0, we must have '^{hN{z) — h{z)) G (— arccos(l — 
e), arccos(l — e)). This implies that hN{z) converges uniformly to h{z) 
on \z\ < r. Since r is an arbitrary positive number less than R, we see 
that (2) holds. 

That (2) implies (1) is an immediate consequence of the equality 
^h„{z) ^ n^^i(l + efc;z^)"^ 

The equivalence of (2) and (3) follows from the equality h'^{z) = 
Ylk=i ^kkakZ^~^ / {1 + tkZ^) and the observation that if the series in (3) 
converges on \z\ < R then it converges compactly. 

Finally, the equivalence of (3) and (4) is a consequence of the in- 
equalities \kakZ^~^\/2 < \ekkakZ^~'' / {I + ekZ^)\ < \kakZ^~''\/ {I - \z\) 
and the remark that if the series in (3) or (4) converges for \z\ < R 
then it converges absolutely and compactly. 

The absolute convergence of the series in (3) and (4) is clear. The 
absolute convergence of the series in (2) follows from the absolute con- 
vergence of the series (4) and the inequality |Log(l + z)\ < 2\z\, which 
holds for all sufficiently small z (for example, \z\ < 1/2 works). □ 
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Lemma 2.3. Suppose that the product f{z) = U'^^^^l + enz'^)°'" con- 
verges compactly on \z\ < R for some R such that < R < 1. Then 
/W(0) = Oforl<i<k. 

Proof. Let h{z) = X]^fc'^nLog(l + e„2;"). By Proposition 12.21 the 
series on the right converges compactly to h, and f = e^. It is clear 
that h^^\0) = for 1 < i < k, as each summand has this property. 
The lemma follows now by differentiation of the equality f = e^. □ 

Proposition 2.4. Suppose that the product f{z) = n^i(l + e„2;")°", 
where ti = ±1, converges compactly on \z\ < R for some R > 0. Then 
f{z) = 1 + Y.7=i ^nZ"-, where 



and the summation ext 
such that ki + 2k2 + . . . + nkn = n. 



and the summation extends over all non-negative integers ki, k2, ■ ■ ■ , kr, 



Proof. By Lemma 12. 3[ the coefficient 6„ coincides with the coefficient 
at z"" in the Taylor expansion of n^=i(l + ekz'')"K The result follows 
now from the binomial series expansion and the Cauchy formula for 
multiplying power series. □ 

Corollary 2.5. Suppose that the sequences bi,b2, ■ ■ ■ and ai, a2, ■ ■ ■ are 

related by ([1]). Then all the numbers 6i, &2; • • • ore integers if and only 
if all the are integers. 

Proof. Since (^) are integers for any integers a and k, the integrality of 
the Oj's implies the integrality of the 6j's. Note that in the formula (1) 
for bn the only contribution of a„ is the monomial e„a„. This observa- 
tion and a straightforward induction on n show that the integrality of 
6j's implies the integrality of a^'s. □ 

Lemma 2.6. Let {Hn) be a sequence defined recursively by Hi = 1 and 
Hn = Hci. Then < Hn < n"^ for every n. In addition, if (bn) 

d\n,d<n 

is any sequence and the sequence (an) is defined by an = J2d\n ^dHn/d, 
then bn = an- ^ a^. 

d\n,d<n 

Proof. Define a sequence (J„) by Ji = 1 and J„ = — 1 for n > 1. The 
definition of (Hn) is equivalent to the equality (Hn) * (Jn) = (En), 
where * denotes the Dirichlet convolution, Ei = 1, and En = for 
n > 1. This means that (Hn) and (J^) are inverses of each other 
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under the operation *. Thus, if (a„) = * (Hn), then (6^) = (a„) * 
(J„). This proves the second part of the lemma. For the first part, 
recall that Ed|n^' < ^'Ilp\ni^ - P~')~^ for s > 0. It follows that 
Ylid\n — ""-^Cls), where ( is the Riemann zeta function. Let s be such 
that ({s) < 2. Then we claim that Hn < n^. Indeed, this is clear 
for n = 1. Assuming that it holds for indexes less than n, we have 

Hn = Ed\n,d<nHd < Ed|n,d<n < C(s)«' - < . ThuS OUr ckim 

follows by induction. Taking s = 2 yields the first part of the lemma, 
because C(2) = n'^/G < 2. □ 

Lemma 2.7. Let g{z) = Yl'^=i dn^"^'^ be analytic in the circle \z\ < R 
for some < R < 1. Let Cj = ±1. Then there exists a unique sequence 
(ofc) such that 

oo 

(2) g{z) = J2 enna^z^-'/il + e^z^ 

n=l 

for all z such that \z\ < R. Moreover, 

(3) -gn = Y,dc^d{-edr^''-'. 

d\n 



Proof. Note that 



1 + enz^' 



k=l 



Suppose first that ([2]) holds. Then the right hand side of (Ej) converges 
compactly on \z\ < R and therefore 



oo oo 



Y,9nZ- = -zg{z) = Y,Y.^ 

n=l n=l k=l 



^^dadi-ed)-'-^ ^z" 



n=l d\n 

Comparing the coefficients at we get the formulas ([3]). A straight- 
forward induction establishes that for any sequence there is a 
unique sequence (q;„) such that holds for all n. This shows the 
uniqueness. For the existence, it remains to show that the series 
^1^^ e„na„2;"~V(l + ^n^") converges for the sequence (q;„) defined 
by ([3]) and all z such that 1^1 < R. By Proposition 12. 2[ it suf- 
fices to show that ^a^-^" converges for \z\ < R. We claim that 
\nan\ < ^d\n\9<i\Hn/d, whcrc (Hn) is the sequence from Lemma [2l6l 
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Indeed, for n = 1 this is clear. Assuming that it holds for indexes less 
than n, we see that 

\nan\ < \gn\ + ^ d\ad\ < \gn\ + ^ ^ \9e\Hd/e = ^ \gd\Hn/d 
d\n,d<n d\n,d<n e\d d\n 

(the last equality holds by Lemma 12. 6p . Thus our claim follows by 
induction. Since Hn < by Lemma [2.61 we see that 



(4) \cyn\<nJ2 



9k\ 



k=l 



Since R < 1 and J2n9n^"' converges for \z\ < R, each of the follow- 
ing series also converges for \z\ < R hj standard properties of power 
series: En l^?n|^^ En(ELi En ^(ELi l^?^!)^"- By©, the 

convergence of the last series implies the convergence of En <^nZ^- D 

Remark 2.8. The convergence of En '^n^"' can be obtained in a differ- 
ent way as follows. Let be the sequence defined by —gn = Ed[n 
(which is the sequence (a„) obtained when en = — 1 for all n). Then 
the inequality |nd;„| < Ed|n \9d.\ follows easily from Mobius inversion 
formula. This implies the convergence of En*^"^"- Now note the fol- 
lowing identity: 



1 - z*^ 1 + e„;z" 2 I- z 



2n 



Using this formula, we can rewrite the series E^i —nanZ'^^ / (1 
term by term, starting with n = 1, into E^i ^nnanZ^~^ / (1 + e„2"). It 
is not hard to see that for n = where m is odd, the a„ obtained 
in this way is of the form ±0^1 ± . . . ± Amt, where rrii = 2*^m and 
si < S2 < ■ ■ ■ < St < s. This observation and the convergence of 
En«„2;" easily imply the convergence of En'^"-^"- 

Theorem 2.9. Let f{z) be an analytic function without zeros in the 
disk 1^1 < i? < 1 and let (e„) be a sequence of ±1. Then there exists 
a unique sequence (a„) such that the product /(O)nj^^(l + e„2")"" 
converges compactly to f on \z\ < R. Moreover, if f{z) = /(0)(1 + 
E^i ^n^") f'{z)/ f{z) = Yl'^=i 9nZ"'^^ then the following formulas 
hold: 



(5) &n = 5^ 



^1 ^2 • • • ^n ' 
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n-1 

(6) nbn = gn + ^ hgn-k, 

k=l 

(7) -^?„ = $^rfad(-Q)"/''^^ 

d\n 

Proof. We may assume that /(O) = 1. The function g{z) = f'{z)/f{z) 
is analytic in \z\ < R. By Lemma [2.71 there exists a unique sequence 
(ctfe) such that 

oo 

(8) g{z) = J2 enria„z"-V(l + e„^") 

n=l 

for all z such that \z\ < R. By Proposition \2.2\ we get log/(z) = 
ZZi «nLog(l + en^") and /(^) = B^ ^(l + e„z")"". 

Formula ([S]) has been obtained in Proposition I2.4[ Formula (E]) fol- 
lows from the equahty /' = fg, i.e. from 

oo / oo \ oo 

n=l \ n=l / n=l 

Finally, ([7]) has been established in ([3]) of Lemma 12.71 □ 

Remark 2.10. The three natural choices for the sequence (e„) are 
e„ = — 1 for all n, e„ = 1 for all n, and e„ such that a„ has non- negative 
real parts for all n. That the third choice always exists follows easily 
from the rewriting procedure described in Remark 12. 8[ Unless some of 
the cXnS are purely imaginary, such a sequence (e„) is unique. 

Remark 2.11. Let (e„) be a sequence of ±1. Starting with a function 
f{z) = 1 + J2'^=i^nZ"', analytic and without zeros in \z\ < R, we can 
compute the exponents a„ recursively in any one of the following ways: 

(i) Using formulas ([5]). 

(ii) By Proposition 12. 4[ the sequence is obtained recursively by 
the following rule: en+idn+i is the coefficient at z^~^^ in the Taylor 
expansion of f{z)U'^^^{l + e^z^)'"'^. 

(iii) Using ([7D and ©. 

It is intriguing that Corollary 12 . 51 does not seem to be easily derivable 
just from ([7]) and (P), even though it is a straightforward consequence 
of©. 
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Example 2.12. We apply Theorem 12.91 to the exponential function 
f{z) = and e„ = —1 for all n. Since /'/ / = 1, we see that gi = I and 
Qn = for n > 1. By ([7]) and the Mobius inversion formula we easily 
get a„ = ~jj,{n)/n. Thus we have the following product expansion: 

oo 

(9) e-^ = n(l-^")"^ 

n=l 

which converges for \z\ < 1. Formula (j9]) is not new. It has been stated 
already in [13] (see formula (13) therein). Now it is well known that 
^^^/i(n)/n = (this equality, conjectured by Euler and proved by 
von Mangoldt, is equivalent to the prime number theorem). Thus we 
may write 

(10) »-'=n(T" 

n=l ^ 

Taking z = 1 leads to the equality 



— z 



oo 

e " = I I 77, " 



n=l 



or, equivalently, 

(11) _ ^ ^ ^/i(n)lnn 



1 ^ 

n=l 



Of course, what we did above is just a heuristic argument, as ( fTOl) 
is valid only for \z\ < 1. Nevertheless, ( ITTi) is correct and it has been 
stated by Mobius [13j (who used heuristic arguments similar to ours, see 
his formula (21)) and proved by E. Landau [9j. We hope that the above 
heuristic argument provides evidence that the product decomposition 
established in Theorem 12.91 may be a source of interesting results in 
number theory. Yet another application will be discussed in the next 
section. 

Example 2.13. We apply Theorem 12.91 to the exponential function 
f{z) = e~ and e„ = —1 for all n. Since /'// = —1/(1 — = 
Y^'^=i{~^)^"'~^ y see that gn = —n for all n. By ([7]) and the Mobius 
inversion formula we easily get a„ = (f){n)/n. Thus we have the follow- 
ing product expansion: 



(12) e^ = n(l-^") 

which converges for \z\ < 1. 



0{n) 

i-z-^) " 

n=l 
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3. Arnold's Conjegture 

The recursive formulas ([6]) are often called formulas of Newton. More 
precisely, substituting g„ = —gn we get 

(13) Qn + biq-a-i + . . . + bn-iQi + ub^ = 0. 

Newton observed in his Arithmetica Universalis, published in 1707, 
that when —bi, . . . , (—1)^6^ are the elementary symmetric functions of 
Xi, . . . ,Xk (and 6„ = for > k) and g„ = a;" + . . . + then the 
relations f[T^ hold. Perhaps a bit less known are the following explicit 
formulas, which (in the case of symmetric polynomials) go back to 
Girard (1629) and Waring (1762): 



U4j qn-n2^[ i) kr\k2\...kn\ ' ' " 



(-l)^'i+^2+...+fe„ fc fc, fc„ 

= ^ k,\k,\...kj It; \-2) •••In J • 

where, in both formulas, the summation extends over all non-negative 
integers fci, k2, . . . ,kn such that ki + 2^2 + . . . + nkn = n. See [7j for 
more about these formulas. As observed by Moree in [12] (and by many 
others before), Newton's formulas relating the symmetric functions and 
the power sums follow easily from (E]). In fact, if f{z) = 1 + biz + . . . + 
bkz'^ is a polynomial, then f{z) = {1 — Xiz){l — X2z) ... (1 —Xkz), where 
X\, . . . , x/j; are the roots of the reciprocal polynomial z'^+biz^ ^ + . . .+bk- 
It follows that 



^ ^ j=l ■> n=l 



Thus Qn = — (x" + ^2 + . . . + x^) and Newton's result follows from ([6]). 

Now let us apply (I7j) with e„ = — 1 for all n. By the Mobius inversion 
formula, we get 



(16) nan = -y^Jdfijn/d). 

d\n 



Assume now that bi, . . . ,bk are integers. Then, by Corollary 12. 5[ all a„ 
are integers too. Therefore we get the following result: 
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Theorem 3.1. Let xi, . . . ,Xk be the roots of a monic integral polyno- 
mial q{x) of degree k and let g„ = x" + X2 + . . . + x^. Then 

^^grfyu(n/(i) = (mod n) . 

d\n 

Applying Theorem 13 . 1 1 when q{x) is the characteristic polynomial of 
an integral k x k matrix A and n = is a power of a prime p, we get 
the following result. 

Theorem 3.2. Let A be an integral k x k matrix and n = be a 
power of a prime p. Then trA^"^ = trA^"^ (mod p™) . 

Theorem 13.21 has been conjectured by Arnold ([2], [3]), and it has 
been the subject of several recent publications ([H], [18], [19], [TT]). 
even though it can be found in papers going back to the 1920s ([8], 
|16]). In [11] we proved a more general result using a different method. 
However, the methods developed in the present paper lead naturally to 
a discovery of Arnold's conjecture and the resulting proof is short and 
aesthetically pleasing. We should mention that already the paper of 
Moree [T2] contains a similar proof of f ll6p . though the above arithmetic 
consequences of this equality have not been addressed there. 

4. Numerical Method 

In this section we describe the numerical method mentioned in the 
introduction. We denote by pn the nth prime number. Let / be a 
function analytic and non-zero in the closed disk 1^1 < -R < 1, /(O) = 1, 
/'(O) = 0. Let m be such that Rpm > 1- Our goal is to approximate the 
product nfclm/(l/Pfc) to high accuracy, as many constants in number 
theory appear in such a form. Our strategy here is very similar to the 
one developed by Moree [12] in the special case when / is a rational 
function satisfying some additional properties. 

By Theorem 12. 9[ there is a product decomposition 

00 

f{z) = i[{i-z-r". 

n=2 

Thus 

00 00 00 00 00 00 

n /(i/pfc) = n 11(1 -p^"")" = n wi^-p-t'^r- = nc-H^"- 

k=m k=m n=2 n=2 k=m n=2 

where Cm(s) = nr=m(l -Pfc ')"^ = C(s) nr=/(l-Pfc') is a partial zeta- 
function (note that the change in the order of multiplication is allowed 
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as the product is absolutely convergent). The key observation behind 

oo 

our method is that the product ]^Cm(^)~°" converges rapidly. More 

n=2 

precisely, we have the following result. 



Theorem 4.1. Let f{z) 



n 



n=2\ 



2;")"" compactly converge in 



the disk \z\ < R < 1. Let B be an upper bound for \f'(z)/f{z)\ on 
\z\ = R. Let m be such that Rpm > 1 d'^d let M > m be such that 
[e - l)Bpm 1 



C{R,B,m,M) :-- 



Rpm - 1 {Rpn 



< 1. Then 



17) 



M 



[n 



k=m 



n=2 



< C{R,B,m,M) 



M 



n=2 



(n 



Proof. Let /'/ / = J2'^=i 9nZ"'~^ and. By Cauchy's inequality, we have 
Ifi'nl < B/R^ for all n G N. By formula ([7]) and the Mobius inversion 
formula, we have nan = — J^dinddf^i^/d)- It follows that 



I ^\ - J^n 



In addition, 



00 



1-n 



k=p„ 



for all n > 3. Using these estimates and the inequality ln(l + x) < x 
(for a; > 0) we get 



an In Cm (n) 



n=M+l 



B 



< 



Pn 



n=M+l 

Bp„ 



1 



n=M+l (^P"^)" ^P"^ ~ ^ {RPn.y 

Note now that for 1 2:1 < 1 we have |1 — e^l < (e — It follows that 



n 

n=A/+l 



in 



< 



l)Bpr, 



1 - exp I E Crr 

n=M+l 

C{R,B,m,M) 



in 



< 



Rpm-l {RpmV 

provided M is such that C{R,B,m,M) < 1. Since | UT=m fC^/Pk)] = 
I n^^2 Cm(^)~""|! ths theorem follows now by multiplying the last in- 



equality by I Un=2Cniin) ""|. 



□ 
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The following two examples illustrate our method. 

Example 4.2. We will compute the first fifty decimal digits of the 
following constant from the paper of Ramanujan [15]: 



Ai = lim 



Inn 



n— )-oo n 



" 1 -j^ oo / 



where d{k) is the number of divisors of k 
-\n{l-z 



Let f{z) 



VT^. Then v^Ai = nr=i/(VPfc)- The 



function / is holomorphic and non-zero in the unit disk. A straightfor- 
ward computation yields 

f (^) _ -1 1 1 

f{z) (l-^)ln(l-^) z 2{l-z)- 

Using the inequality \{1 — z) ln(l — -z)] > — {z]"^ we get 



fiz) 



< 



1 

+ ^ + 



2(1 



2 3 



\z\) \z\ 2(1-1^1)- 

This gives an estimate \ f'/f \ < 18 for \z\ = 0.9. Thus we may take 
R = 0.9, B = 18, and m = 7, so = 17 and Rpm > 15. Take M = 48. 
Then 

C{R,B,m,M) < ^ < 10"^^ 



Thus, by Theorem 14.11 



15^8 



, 6 48 

^1 - 7i n f(^/p^^ n 

k=l n=2 



-54 



6 48 

/^n/(vp.)ncrH 

k=l n=2 



< 10 



Now we calculate the exponents a„, n = 2, 3, . . . , 48 and the product 
11^=1 /(l/Pfc) nn=2C7(^)~"" accurate to 54 decimal digits and get 
Ai = 0.54685595528047446684551710099076178991021048592974 . . . (the 
computation has been done with 211 accurate digits using PARI/GP). 

Example 4.3. Let d{k) and denote the number and the sum of 
divisors of k, respectively. Then a{k)/d{k) is the average divisor of k. 

1 " 

The average of the average divisor is then the quantity — } a(k)/d(k). 



k=l 



In [5l Thm. 4.1] it is proved that this quantity is asymptotically equal 
to cn/ Inn, where 



^ oo 



3/2 
Vk 



\ Vpk 



In 1 



Pk 



^ oo 

7^11 /(VP.). 

k=l 
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where f{z) 



ln(l + z) 



V. I. Arnold, in his recent book attributes 



Zy/1 — Z 

this asymptotic to A. A. Karatsuba (note however that the formula for 
Ci = c in the footnote on page 78 of [4J is incorrect: it has a factor 
of l/vr instead of I/a/tt). According to Arnold, M. Korolev computed 
c ~ 0.7138067 .... We will see that only the first 5 digits are accurate. 
The function / is holomorphic and non-zero in the unit disk and 

f (^) _ 1 1 , 1 

fiz) {l + z)\n{l + z) z 2{l-z)' 

The same estimates as in Example 14.21 allow us to take R = 0.9, B = 18, 
m = 7, and M = 48 and get 



6 



48 



k=l n=2 

Now we calculate the exponents a„ 
r6 x/i / \ tt48 



< 10 



-54 



6 



48 



TT 



fe=l n=2 

n = 2, 3, . . . , 48 and the product 
nfc=i nn=2C7(^)~"" accurate to 54 digits and get 

c = 0.71380993049991415224401060402799291827213336525147. . .. 
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